Thermal Relaxation of Charm in Hadronic Matter 
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The thermal relaxation rate of open-charm (D) mesons in hot and dense hadronic matter is 
calculated using empirical elastic scattering amplitudes. D-meson interactions with thermal pions 
are approximated by D* resonances, while scattering off other hadrons (K, 77, p, uj, K* , N, A) 
is evaluated using vacuum scattering amplitudes as available in the literature based on effective 
Lagrangians and constrained by realistic spectroscopy. The thermal relaxation time of D-mesons 
in a hot n gas is found to be around 25-50 fm/c for temperatures T=150-180 MeV, which reduces 
to 10-25 fm/c in a hadron-resonance gas. The latter values, argued to be conservative estimates, 
imply significant modifications of D-meson spectra in heavy-ion collisions. Close to the critical 
temperature (T c ), the spatial diffusion coefficient (D s ) is surprisingly similar to recent calculations 
for charm quarks in the Quark-Gluon Plasma using non-perturbative T-matrix interactions. This 
suggests a possibly continuous minimum structure of D 3 around T c . 

PACS numbers: 12.38.Mh, 24.85. +p, 25.75.Nq 

Keywords: Charm transport, hot hadronic matter, ultrarclativistic heavy-ion collisions 
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I. INTRODUCTION 

After the initial discovery of a new state of matter in 
high-energy nuclear collisions at the Rclativistic Heavy 
Ion Collider (RHIC), the focus is now shifting to quan- 
tifying the properties of what is believed to constitute 
a strongly coupled Quark-Gluon Plasma (sQGP). Ba- 
sic quantities characterizing the medium are its thermal 
spectral functions and transport properties. In heavy-ion 
collisions (HICs), the former are most directly studied 
in the electromagnetic (vector) channel via the thermal 
emission of lepton pairs (cf. Ref. [l| for a recent review). 
The latter, however, are best studied using observables 
with small but controlled deviations from thermal equi- 
librium. Thus, heavy quarks (charm (c) and bottom (6)), 
whose thermal equilibration time is expected to be of the 
order of the QGP lifetime in HICs, are a promising tool to 
quantify flavor transport, and eventually deduce general 
properties of the QGP as formed in these reactions Q. 

The large masses of c and b quarks (m Ct b) enable us 
to assess the modifications of their momentum spectra 
in HICs via a diffusion process in an evolving back- 
ground medium as formulated, e.g., within a Fokker- 
Planck equation Q (typical early temperatures of the 
medium produced at RHIC, T ~ 250 MeV 0], are well 
below m c ,b — 1.5,4.5GeV). A reliable determination of 
the heavy-quark (HQ) transport coefficients in the QGP 
depends on several components. First and foremost these 
are microscopic calculations of the thermal HQ relaxation 
rate in the QGP [IlHlJ] ( see Ref - II for a review). Sec- 
ond, the coefficients need to be implemented into a re- 
alistic bulk medium evolution (see, e.g., Ref. [l5| for a 
recent discussion). Third, heavy-flavor (HF) interactions 
in evolution phases other than the QGP have to be evalu- 
ated, i.e., in the so-called prc-equilibrium phase as well as 
in the hadronic phase. The former is of a relatively short 



duration, Ar pro < 1 fm/c, and is sometimes mimicked by 
reducing the formation time of the QGP. The duration 
of the hadronic phase is substantially longer, Arhad— 5- 
10 fm/c. Its relevance for HF phenomenology is further 
augmented by the fact that the hadronic medium inherits 
the full momentum anisotropy from the QGP, believed to 
be close to the finally observed one. Thus, even a rather 
weak coupling of HF hadrons to hadronic matter can lead 
to noticable contributions to their elliptic flow. Further- 
more, if the QGP realizes a minimum in its viscosity-to- 
entropy-density ratio, rj/s, close to the (pseudo-) critical 
temperature, T c ~ 170 MeV, a hadronic liquid close to 
T c should possess similar properties. This is usually re- 
ferred to as a "quark-hadron duality", as suggested, e.g., 
in calculations of thermal dilcpton emission rates [l| . 

Charm diffusion in hadronic matter has received little 
attention to date (see Ref. [l6| for a very recent estimate 
using heavy-meson chiral perturbation theory). Its po- 
tential relevance has been noted in Ref. 0] based on cal- 
culations of D-meson spectral functions in nuclear matter 
using effective hadron Lagrangians [l7], EH , as well as in 
a hot pion gas fX9| . In the present paper, we augment 
these works to evaluate charm diffusion in hadronic mat- 
ter. Since the latter features many resonances at tem- 
peratures approaching T c , we not only utilize Dit and 
DN interactions but also scattering amplitudes off ex- 
cited hadrons (K, r\, p, uj, K* , A), as constructed in the 
literature using effective Lagrangians and constrained by 
charm-resonance spectroscopy. In this sense we provide 
a lower estimate of the diffusion coefficient, based on ex- 
isting elastic D-hadron amplitudes. 

Our paper is organized as follows. In Sec. [XT] we "re- 
construct" microscopic models for Dtt scattering in a hot 
pion gas (via s-channel resonances; Sec. Ill A| ). for D scat- 
tering off strange and vector mesons (Sec. Ill Bj) , and off 
baryons (Sec. Ill Cj) . by parameterizing pertinent scatter- 
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ing amplitudes. In Sec. IIIII these are applied to calcu- 
late thermal Z?-meson relaxation rates and diffusion co- 
efficients in hot hadronic matter at vanishing chemical 
potential, first in a pion gas (Sec. IIII A|) and then in a 
resonance gas (Sec. IIII B[) . and finally including chemi- 
cal potentials as appropriate for heavy-ion collisions at 
RHIC (Sec. IIII C|) . Conclusions are given in SecHVl 



II. 



D-MESON SCATTERING AND WIDTH IN 
HOT HADRONIC MATTER 



In this section we recapitulate basic elements of the 
_D-hadron scattering amplitudes and apply them to cal- 
culate pertinent D-meson widths in hadronic matter. We 
only employ known amplitudes from the literature for the 
most abundant hadrons in a hot medium and combine 
them into a lower-limit estimate for the D-meson width 
(and D- meson relaxation rate in Sec. IIII[) . We first focus 
on pion-gas effects, followed by interactions with strange 
and vector mesons, as well as hot nuclear matter. 



A. Hot Pion Gas 

Following Ref. D interactions in the pion gas are 
dominated by its chiral partner, Dq(2308), by the vec- 
tor D*(2010) and by the tensor meson 1^(2460). Phe- 
nomenological control over their interaction vertices has 
become possible due to new observations of £>-meson res- 
onances by the BELLE Collaboration [2(|. Especially, 
the large width of the DS(2308), Fjf = 276 ± 99MeV, 
attributed to 5-wave pion decay, leads to a large DttDq 
coupling constant. With the constituent-quark model as- 
signment of isospin 1—1/2 for £>-states, the pertinent for- 
ward scattering amplitudes have been parameterized in 
Breit-Wigner form as 



M(s,e = o) 



E 

3=0 



8n^ (2j + l)VSP 



s M] 



lot 

3 



(1) 



where y/s and k denote the center-of-mass energy and 3- 
momentum in the Dtt collision, respectively, and j is the 
resonance spin. The total resonance decay width T tot , is 
assumed to be saturated by the partial one, 7. for D* 
and Dq, while for the narrow state D\ (r| ot = 45.6 ± 
12.5 MeV) a branching fraction estimated by the particle 
data group is employed [2l[ . The resulting total Dir cross 
section is displayed in the upper panel of Fig. [TJ 

The D-meson self-energy in a pion gas can now be 
obtained by the standard procedure of closing the in- 
and outgoing pion lines of the forward Dtt amplitude 
with a thermal pion propagator. In the narrow- width 
approximation for the pion propagator, the D-meson self- 
energy takes the form 



n(p D ;T) 
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FIG. 1: (a) Total elastic cross sections for D-tt (solid line), S- 
wave D-K (dashed line) and D-p (dash-dotted line) scattering 
based on the invariant amplitudes as discussed in the text; 
vertical lines indicate the respective thresholds, (b) Total D- 
meson scattering width in a meson gas (dash double dotted 
line), consisting of contributions from Goldstone bosons and 
vector mesons. 



where fs is the thermal Bose function and E n the on- 
shell pion energy. Alternatively, one can obtain the colli- 
sion rate (or on-shell width), T=—IiJiIi(p 2 D =m 2 D , T)/rriD 
from the Boltzmann equation as 
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where 7^=3 denotes the spin-isospin degeneracy factor 
of the in-medium particle (pion), and p, q and p' \q' are 
the momenta of in- and outgoing particles, respectively. 
Equations ([2]) and (|3|) are related via the optical theorem 
(we have checked their consistency). In the following 
we employ the latter since it is close to the form of the 
thermal relaxation rate discussed in Sec. IIIII below. 

The scattering width of a Z?-mcson at rest in a pion 
gas is displayed in the lower panel of Fig. [T] For T=150- 
180 MeV we find r£,=20-40McV, where the chiral part- 
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ner of the D provides the largest contribution through S- 
wave Dtt scattering. For constant resonance widths, we 
find close agreement with the results plotted in Fig. 2 of 
Ref. , while for energy- momentum dependent widths 
(as quoted in Ref. [HI) our results displayed in Fig. Q] 
turn out to be slightly smaller (by ca. 20%). 

Chiral effective Lagrangians have been applied to D- 
mcson interactions with Goldstonc bosons in Rcfs. (22T - 
[2~H - Once the parameters are constrained by empirical 
information, the resulting scattering amplitudes are very 
similar and closely agree with the resonance ansatz of 
Ref. as employed here. Chiral Lagrangians also pre- 
dict non-resonant interactions. In the (repulsive) S-wave 
1=3/2 Dtt channel the scattering length has been com- 

puted as a D7r =— /(sthr)=— O.lfm in unitarizcd chiral 

1 1 7—3/2 

perturbation theory (%PT) [25[ vs. a D7T =— 0.15fm in 
%PT to NNLO [26|. We employ the former to construct 
the 1=3/2 scattering amplitude as Ai(s) = 8tv y/sf(s). 



where gov is the dimensionful coupling constant and sr 
the complex resonance-pole position. We include the 
three 1=1/2 resonance couplings to Dp and Dui from 
Tab. 7 in Ref. [27j |, two 1=0 and two 1=1 resonances 
with DK* coupling (Tabs. 5 and 4 in (27}, respectively) 
and one 1=0 state coupling to DK* (Tab. 8 in (27) ). As 
a representative, the isospin 1=1/2 Dp cross section is 
shown in the upper panel of Fig. [TJ 

The lower panel of Fig.[T]shows the temperature depen- 
dence of mesonic contributions to the D-meson scatter- 
ing width as calculated from the above amplitudes. The 
width from anti-/kaons is the next largest contribution 
after the pion. The effect of vector mesons is smaller but 
significant, especially for the K*. The total D width in 
a hot meson gas reaches ~80MeV around T~180McV, 
which should be a lower limit since several channels have 
still been neglected, e.g., higher partial waves (except for 
pions) and inelastic channels (e.g. DK* -H> D s ir). 



B. Strange and Vector Meons 



C. Hot Nuclear Matter 



In a hot meson gas, the next abundant species after 
the pions are the strangeness carrying Goldstonc bosons 
and the light and strange vector mesons. 

For D(K, K) interactions wc directly employ the re- 
sults of the (unitarized) chiral effective theory of Ref. 22 1 
(by parameterizing the amplitudes in Figs. 1 and 2 
in there). In the S'-wave DK channel, with isospin- 
strangeness (I, S)=(0,1), the results are constrained by 
the (loosely bound) £>* (2317) state (the analogue of 
_Dn(2308) in Dtt), and again closely agree with Refs. [H, 
[Hjj. The application to the S'-wave DK (I, S) = (I,1) and 
DK (I, S)=(0,— 1) produces Feshbach-type resonances 
(tetra-quarks) right at threshold (£^-=2360 MeV). For 
DK in the (I, S)=(l,— 1) channel the Born amplitude is 
predicted to be repulsive (analogous to 1=3/2 Dtt); the 
scattering length has been calculated as a^~^=— 0.22fm 
in unitarized x?T [25[ vs. a^g=— 0.33fm in in xPT at 
NNLO [Hj]. Analogous to 7=3/2 Dtt, we employ the 
former to construct the pertinent scattering amplitude. 

For Drj scattering we also adopt Ref. [22||, where the 
S-wave (I, S)=(l/2,0) amplitude is governed by a narrow 
state of mass 2413 MeV just below threshold. 

The evaluation of DV scattering (V=p, uj, K*) re- 
quires to go beyond the chiral Lagrangian. This has 
been done in Ref. [2lJ starting from SU(4) flavor symme- 
try and then implementing chiral breaking terms. This 
framework, properly unitarized, recovers the resonance 
poles computed with the chiral Lagrangians, but ex- 
tends to axialvector states coupling to S-wave DV in- 
teractions (in particular L»i(2420), L>i(2427), L> s i(2460), 
-D s i(2536)). A convenient Breit-Wigner parametrization 
of the elastic coupling of DV to the dynamically gener- 
ated resonances has been quoted as 
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To evaluate D scattering off baryons wc follow the 
same strategy as for mesons, employing vacuum scatter- 
ing amplitudes. More elaborate many-body calculations 
for D-mesons in nuclear matter are available [H, [28| , but 
our procedure keeps consistency with the mesonic sector 
and allows for an estimate of the systematic error due to 
in- medium effects (e.g., sclfconsistcncy of selfcnergy and 
in-medium scattering amplitudes). 

We start from Ref. [l8| where the T-matrix results of 
an effective DN interaction with coupled channels have 
been parameterized in analogy to Eq. (0]) as 



\fJ\ 



y/s - m,R + iT R /2 



(5) 



(4) 



where is the resonance mass and Tr the width. The 
T-matrix, Eq. ([5]), is related to the invariant scattering 
amplitude through M(s) = N(s)T(s) with N(s)=(s + 
m 2 N — m 2 D + 2mAr v / i)/(2y / s). The key states are the 
dynamically generated 1=0 A c (2595) and 1=1 E c (2620) 
S-wave bound states. The former is experimentally well 
established, while the latter is ca. 180 MeV too deep com- 
pared to the empirical S c (2800) state. However, the 1=1 
DN cross section above threshold (shown in the upper 
panel of Fig. [2]) is quite comparable to the results of a 
recent meson-exchange model calculation [2^] in which 
the E c is generated at its experimental mass (in fact, at 
threshold the 1=1 scattering length in Ref. [29[ is signifi- 
cantly larger than for the S c (2620) state of Ref. Il3). In 
the 1=0 channel, the scattering lengths of Refs. [18( and 
[29l | agree well. The corresponding cross section is also 
shown in the upper panel of Fig. [2] 

The DN scattering amplitude can be inferred from 
DN due to C-symmctry of strong interactions. We have 
found no evidence in the literature for (multi-quark) reso- 
nances in this system and adopt the D~ N elastic S-wave 
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of Fig. [2]). When adding N and A contributions, the 
baryon-induced D-meson width computed here amounts 
to ca. 15 MeV at T~180MeV. Note again that we have 
neglected higher partial waves as well as higher excited 
resonances including strange anti-/baryons. 

It is instructive to compare our nucleon-induced width 
to a sclfconsistent many-body calculation [l?}- From 
Fig. 6 in Ref. [H we read off r r ,~100(80) MeV at 
T=100(150)MeVand g N =Qo, vs. r D ~75(65) MeV in our 
approach. This indicates that neglecting in-medium ef- 
fects does not lead to an overestimate in our calculation. 



III. THERMAL RELAXATION IN HADRONIC 
MATTER 
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FIG. 2: (a) Elastic cross sections for S-wave D-N scattering 
in the 7=1 (solid line) and 1=0 (dotted line) channel, and D- 
A scattering (dashed line, 1=1; vertical line: DA threshold) 
based on the invariant amplitudes as discussed in the text; (b) 
D-meson widths from D-N scattering (solid line, including 
both 1=1 and 1 = channels) and from D-A scattering 
(dashed line) at vanishing chemical potential. 



amplitude calculated in Ref. [18|. The pertinent scatter- 
ing and relaxation rates are about a factor of ^3 smaller 
than from UJV scattering. 

The only available calculation of DA we are aware 
of has been conducted in Ref. [HJ, within the same 
framework as our DN amplitudes are based on. In 
the only available 1=1 S'-wave channel the parametriza- 
tion, Eq. (O, reflects a rather deep bound state 
(?n#=2613MeV, r^O, 5=8.6). The cross section is 
shown in the upper panel of Fig. [2] Unlike the DN case, 
the DA system is predicted to support a shallow 1=1 
bound state (TO fi =2867McV, T R ~0, 5=5.8). As a re- 
sult, the contribution of DA scattering to the D-mcson 
width and relaxation rate is about half of that from DA 
scattering. 

The lower panel of Fig. [2] shows that the width of a 
D-meson at rest from scattering off thermally excited nu- 
clcons and A's at vanishing chemical potential (/Lts=0) is 
comparable to that of light vector mesons (cf. lower panel 



The standard expression for the thermal relaxation 
rate of a particle (D) in a heat bath in terms of its scat- 
tering amplitude on medium particles (h) reads Q 



1 d 3 q d 3 p> dV h 
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with p (q) and p' (<f' ) being the D-meson (h) momentum 
before and after the interaction, respectively. The form 
of this expression is very similar to the one for the total 
width, Eq. ([3]). The latter can be expressed in terms of 
the average defined above as T = ((1)). In Sec. IIII A1 we 
evaluate A(p, T) for D-mesons in a pion gas (h=w) and 
in Sec. IIII Bl for all other hadrons whose amplitudes have 
been constructed in the previous section. In Sec. IIII O we 
evaluate the relaxation rate at finite baryon and meson 
chemical potentials as applicable to the hadronic phase 
below the chemical freeze-out temperature in heavy-ion 
collisions at RHIC. 



A. Pion Gas 

The thermal relaxation rates for a D-meson at rest 
due to scattering off pions in a thermal gas in chemical 
equilibrium (^=0), using the the scattering amplitude 
of Eq. ([1}, is displayed in Fig. [3] as a function of tempera- 
ture. From T=100 MeV to 180 MeV the rate increases by 
about a factor of 7, basically following the increase in pion 
density from 0.2 to lAg . Its magnitude at T=180MeV, 
A ~ 1/25 fm, is small but not negligible. When replacing 
the D7T amplitude in Eq. ([3]) by one yielding a constant 
S'-wave cross section of CT£ 7r =7-10mb, the pertinent band 
for the relaxation rate essentially covers the result of our 
microscopic calculations. The latter are closer to the up- 
per end of the band at lower T but to the lower end 
at higher T > 150 MeV. This reflects the increased ther- 
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FIG. 3: Thermal relaxation rate of D-mesons at momentum 
Pd=0.1 GeV in a gas of pions (solid line), anti-/kaons (dash- 
double-dotted line), nucleons (dash-dotted line) and antinu- 
cleons (dotted line), as a function of temperature at vanish- 
ing chemical potentials using the elastic scattering amplitudes 
constructed in Sec. HU The bands represent results obtained 
with constant 5-wave Dn, DK and DN cross sections of 7- 
10 mb and 10-15 mb, respectively. 



FIG. 4: Thermal relaxation rate of po=0.1 GeV D-mesons us- 
ing empirical amplitudes in a hadron gas at vanishing (solid 
line) and finite (dash-dotted line with squares) chemical po- 
tentials, as well as in a pion gas (dotted line). The dashed line 
corresponds to isotropic D-meson cross sections with mesons 
(7mb) and baryons (10 mb). The filled box at the upper right 
indicates charm-quark relaxation rates in a QGP at 1.2 T c 
from in- medium T-matrix calculations uM. 



mal motion of pions probing higher s/s in the amplitudes 
where the latter decrease. 

Compared to the recent work of Ref. [l6|, where D- 
mcson diffusion in a lukewarm pion gas has been eval- 
uated in heavy-meson chiral perturbation theory, our 
value for the friction coefficient is much smaller; e.g., 
at T=50(100)MeV, Ref. finds A = n/2m D T ~ 
0.00055(0.05)/fm, about a factor of -4(10) larger than 
our pion-gas results, A ~ 0.00015(0.005)/fm. 



B. Hadron Resonance Gas 

A hot hadron gas in equilibrium is characterized by 
an increasing abundance of resonances with rising tem- 
perature. For example, at T=180MeV, the density of 
baryons plus antibaryons is above go and that of mesons 
with masses below 2 GeV is above 3qq. To improve the 
estimate of D-meson diffusion in a pion gas for a more 
realistic hadron-resonance gas, we include rescattering on 
all particles which at T=180MeV and ph=0 have a den- 
sity at least O.lgo, i.e., tt, K, r], p, ui and K*(892) in the 
meson sector and anti-/nuclcons, and A(1232), A(1232) 
in the anti-/baryon sector, using all scattering ampli- 
tudes of Sec. HU The resulting D-meson friction coeffi- 
cient in hadronic matter at vanishing chemical potentials 
increases substantially over the pion-gas result, by about 
a factor of -2(3) at T=150(180) MeV, seeFig.gl The in- 
dividual contributions of K + K and N, N are compared 
to constant-cross-section calculations in Fig. [31 indicating 
that a "constituent" light-quark cross section of 3-4 mb is 



compatible with our lower-limit estimates. A quantita- 
tive decomposition of the individual hadron contributions 
to kinetic D-meson relaxation at T= 180 MeV is given in 
Tab. |H Anti-/kaons provide the next-to- leading contri- 
butions after the pions, while vector mesons, nucleons 
and Deltas play a smaller but non-negligible role. 

As an estimate of medium effects in our vacuum 
Dh amplitudes we have performed a calculation for A 
where we have introduced an in-medium broadening of 
rg od =200McV into the Breit-Wigner parameterizations. 
The hnal result for A changes by less than 5%. 

Quantitatively, a value of A~0.1/fm translates into a 
thermal relaxation time of rjj—lOfm. This time scale 
is quite comparable to lifetimes of the hadronic phase 
in Au-Au collisions at RHIC. It is also similar to non- 
perturbative calculations of charm-quark relaxation in 
the QGP using in- medium T-matrix interactions [TJ E3 ■ 
This is encouraging both from a conceptual (in the sense 



Hadrons 


Ll,2J 


A [fm" 1 ] 


7T 


Sl/2,0, -Pl/2,2, Dx/2.4, S'3/2,0 


0.0371 


K + j) 


So,0, Sl,0 


0.0236 


P + uj + K* 


Sl/2,2, 5*0,2, 5*1,2 


0.0129 


N + N 


5o,i, 5i,i 


0.0128 


A + A 


5*1,3 


0.0144 



TABLE I: Contributions to the thermal D-meson thermal re- 
laxation rate at T= 180 MeV indicating the quantum numbers 
of the included scattering channels with L: partial wave, 7: 
isospin and J: total angular momentum. 
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IV. SUMMARY AND CONCLUSION 
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FIG. 5: Spatial D-meson diffusion coefficient in units of the 
medium's thermal wavelength, l/(2nT). The filled box at the 
lower right indicates the range of values calculated for charm 
quarks in the QGP at 1.2 T c within an in-medium T-matrix 
approach 



of a quark-hadron continuity through T c ) and a prac- 
tical point of view (relaxation rates of this magnitude 
currently provide a fair phe nomenology of current heavy- 
flavor data at RHIC 

In Fig. [5]we display the spatial D-meson diffusion co- 
efficient, D s = T I (mi)A(j) = 0,T)), in hadronic matter. 
When normalized to the thermal wavelength, l/(27rT), 
this quantity decreases with T, reaching a value of ^5 
at T=180McV. Again, this is surprisingly close to T- 
matrix results for charm quarks in the QGP [lH, and, 
together with those results, suggests a minimum across 
the hadron-to-quark transition. 



C. RHIC Conditions 

In relativistic HICs the chemical freeze-out of hadron 
ratios [3l[ at a temperature of T chem ~170 MeV is signifi- 
cantly earlier than thermal freeze-out of the light hadrons 
at Tf o ~100 MeV. Therefore, to conserve the observed par- 
ticle ratios in the hadronic evolution, effective chemical 
potentials are required, reaching appreciable values even 
at RHIC energies (H, e.g., /i 7r (T=100MeV)~80MeV. 
We implement the chemical potentials into the ther- 
mal hadron distribution functions and recalculate the D- 
meson equilibration rate, Eq. ([5]). As a result, the latter 
is enhanced at temperatures below T c h em , staying above 
1/(25 fm) for T>130MeV (cf. Fig. g]), implying notice- 
able modifications of D-meson spectra in the hadronic 
phase of nuclear collisions at RHIC. For example, if the 
hadronic evolution lasts for Arhad— 5 fm, the expected 
modification amounts to ca. (1 — exp[A Arhad]) — 20%. 



We have evaluated kinetic transport of D-mesons in 
hot hadronic matter by elastic scattering off the 10 most 
abundant hadron species. The interaction strength with 
mesons and baryons has been estimated from existing 
microscopic models for D-hadron scattering, constrained 
by chiral symmetry and vacuum spectroscopy. In a pion 
gas at T=100McV, Dn resonance scattering leads to a 
relaxation rate which is substantially smaller than what 
has been found in a recent calculation using heavy-meson 
chiral perturbation theory. Yet, when extrapolating our 
full results to temperatures in the vicinity of T Cl the re- 
laxation rate reaches ca. 0.1/fm, translating into a spatial 
diffusion coefficient of D s ~ 5/(2irT). This is compara- 
ble to non-perturbative T-matrix calculations of charm- 
quark relaxation in the QGP. On the one hand, this 
suggests a rather smooth evolution of charm transport 
through T c , i.e., a kind of "duality" of hadronic and 
quark-based calculations, reminiscent of what has been 
found for dilcpton emission rates. On the other hand, it 
implies that quantitative calculations of D-mcson spec- 
tra in heavy-ion collisions have to account for hadronic 
diffusion. This insight is reinforced once chemical freeze- 
out is implemented into the evolution of the hadronic 
phase (via effective chemical potentials), with an esti- 
mated modification of D-meson observables of at least 
20%. The apparent agreement of hadron- and quark- 
based approaches, when extrapolated to around T c , is 
encouraging, especially since the magnitude of the trans- 
port coefficient is compatible with the phenomenology of 
current heavy- flavor observables at RHIC. Our findings 
thus pave the way for an improved theoretical accuracy 
which will be needed to take advantage of upcoming pre- 
cision measurements at RHIC and LHC. 

Note added. Two subsequently submitted papers 
have also addressed hadronic D-meson diffusion. In 



Ref. [33[ the use of unitarized chiral effective D7r inter- 
actions leads to relaxation rates in a pion gas in close 
agreement with our results. In Ref. |34j D-hadron in- 
teractions have been evaluated using Born amplitudes, 
leading to relaxation rates significantly larger than our 
results. 
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